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SUMMARY 

An analog investigation of several turbojet control configurations 
was made. Both proportional and proportional-plus -integral controllers 
were studied, and compensating terms for engine interaction were added 
to the control system. Data were taken on the stability limits and the 
transient responses of these various configurations. Analytical expres- 
sions in terms of the component transfer functions were developed for 
the configurations studied and the optimum form for the compensation 
terms was determined. 

It was found that the addition of the integral term, while making 
the system slower and more oscillatory, was desirable in that it made 
the final values of the system parameters independent of source of dis- 
turbance and also eliminated droop in these parameters . 

Definite improvement in system characteristics results from the 
use of proper compensation terms . At comparable gain points the com- 
pensated system is faster and more stable. Complete compensation elimi- 
nates engine interaction, permitting each loop to be developed to an 
optimum point independently. 


INTRODUCTION 

Turbojet engines with a fixed-area exhaust nozzle do not present 
too difficult a control problem because only one input variable, fuel 
flow, is manipulated to maintain desired engine speed or temperature. 

A single closed-loop system, incorporating overspeed and overtempera-, 
ture protection along with a schedule of fuel flow to prevent surge on 
acceleration, will accomplish the necessary control function. When a 
variable-area exhaust nozzle is added to such an engine, however, the 
control problem becomes more complex because two input variables are 
available and these should be so controlled that the engine is at aH 
times operating in a safe and efficient manner. When more than one 
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input variable to an engine is controlled, the resulting system is a 
multiple-loop configuration. A general discussion of multiple-loop 
systems with a specific exainple of an aircraft reciprocating engine con- 
trol is given in reference 1. 

In the specific case of an engine where speed and temperature are 
to be controlled by manipulation of fuel flow and exhaust-nozzle area, 
it follows that two possible double-loop systems can be suggested. In 
one case, speed can be controlled by exhaust -nozzle area, and tempera- 
ture, by fuel flow. In the second system speed can be controlled by 
fuel flow vhile temperature is controlled by exhaust -nozzle area. A 
basic characteristic of turbojet engines is that a change in fuel flow 
or area causes both speed and temperature to change. Therefore, when- 
ever these engine parameters are used in a double-loop control config- 
uration, a disturbance in one loop will introduce an error signal into 
the other loop. This characteristic will be referred to herein as the 
interaction effect which exists between the individual control loops in 
a double-loop system. As a result of such interaction, an unstable 
system, or one having very oscillatory responses in some regions of con- 
trol operation, can result even though each loop may be inherently stable 
when used alone. In order to stabilize a system of this form, it gen- 
erally becomes necessary to reduce the loop gains or sensitivities; but 
this is accomplished at the expense of an increase in response time for 
the conplete system. 

A general algebraic method of analysis has been applied to the de- 
termination of control requirements for multiple -loop engine control 
systems and is presented in reference 2. It was shown therein that con- 
trol systems could be designed so as to be noninteracting, which means 
that each loop in such a configuration can then be considered as acting 
independently in the combined system. Further analysis also indicates 
that a noninteracting control system will offer advantages in regard to 
improved stability and faster response than are possible with the cur- 
rent interacting configurations. It was considered important, there- 
fore, to have an understanding of both interacting and noninteracting 
double-loop systems because the more complex engine types being developed 
at present, along with the demands for faster responding power plants, 
necessitate the use of such systems . 

For this reason an investigation was initiated at the NACA Lewis 
laboratory to determine some of the practical aspects of noninteracting 
systems and to compare these with an interacting configuration. Sta- 
bility limits and response characteristics were obtained for one basic 
double-loop system and also for several modifications of the system. 

An analog computer was used to simulate a current turbojet engine with 
variable jet nozzle along with the necessary sensor and servo compo- 
nents of the engine control. 
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The basic configuration studied is one in vhich speed is controlled 
by fuel flow and temperature by exhaust-nozzle area. One modification 
consisted of adding an integral term to each loop of the system, while 
the second modification consisted of adding a term to compensate in part 
/ for the interaction characteristic of the engine. Stability limits were 

determined for these systems . Three different forms of compensation for 
noninteracting systems were investigated. The investigation was extended 
to present transient response characteristics of the systems to a step 
disturbance in set temperature. The engine is assumed to be operating 
near design speed but at lower than design temperature. An increase in 
thrust is obtained by increasing set temperature. An assumption was also 
made of linearity in the region of the engine operating point. 


COMPUTER AND METHOD 

A high-speed electronic analog computer operating at 4800 times 
real time was used. A number of computational elements of standard fonn 
are available and these can be interconnected by means of plug-in cables . 
A standard square-wave disturbance voltage with a repetition rate of 60 
cycles per second is supplied and, by calibration, its time base repre- 
sents 20 seconds of real engine time. Solutions are presented on a 
group of oscilloscopes so that the transient response of several vari- 
ables due to the applied step disturbance can be observed simultaneously. 
One of the computational elements is a matrix which is used to simulate 
the engine. This method is presented in detail in reference 3. 

A control simulator component is also available which has a trans- 
fer function as follows, where Eq and E^ are output and input 
voltages : 

^ * T^) * ^OS) (5-^) 

where the gain term K, integral, dei-ivative, and lag time constants T^, 
Tq, and T^ are variable. An added feature is that the integral, deriv- 
ative, or lag terms can be switched out if required. The computer also 
contains a number of summing, coefficient, integral, derivative, and lag 
units along with calibration devices which permit a more accurate setting 
of the variables and determination of output voltage values. Provisions 
are also available for photographing the oscilloscope displays . 

A high-speed computer of the type used has the advantage that char- 
acteristic responses over a broad range of possible control settings of 
various systems such as shown on figure 1 can be investigated very 
quickly with minimum effort. (The symbols on fig. 1 and elsewhere are 
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defined in the appendix.) Systems can be quickly changed or modified 
as required by indicated trends of the investigation. 

Stability^ in particular, can easily be determined by the following 
method; With no forcing function or disturbance and with a specific 
value of temperature loop gain set into the computer, the speed loop gain 
can be gradually increased from zero until the entire system becomes un- 
stable, as shown by continuous oscillations of all parameters on the 
oscilloscopes . This procedure can be repeated for a number of values 
of temperature loop gain over the entire range. A plot of the values 
of temperature loop gain versus speed loop gain at which the system be- 
comes unstable can be made from these data, and this curve defines the 
limits of stability for the configuration. When a disturbance is added 
to the system, the transient responses of all pertinent parameters can be 
observed and variations in these responses noted as a function of loop 
gains . 

All engine gain or sensitivity terms used in the simulation were 
nonnalized to rated values and therefore computer output voltages repre- 
senting the transients were proportional to a percent of rated value 
change in all parameters. For the purpose of this report, a 1-percent 
step disturbance was introduced in set temperature. Speed and tempera- 
ture droops (which are defined as the deviation in percent of rated 
value of the parameter in steady state from the desired final value) 
and maximum excursions (which are defined as the maximum deviations in 
percent of rated values of the parameters during a transient, measured 
from the initial starting point) were recorded and plotted as percent 
deviations on the stability limit figures . In addition, the time rises 
(which is defined as the time required to reach maximum excursion) were 
noted and plotted in a similar manner. 

Examination of the resulting maps shows how the transient responses 
vary as a function of both speed and temperature loop gains and also 
permits a rapid comparison of the effect on response that can be ob- 
tained by modifying the system and by using compensation for the inter- 
action normally found in engines. 


SELECTION OF SYSTEMS TO BE INVESTIGATED 

Preliminary analysis of interacting and noninteracting systems was 
conducted to determine the specific configurations to be studied in de- 
tail by analog methods. 

Interacting system . - A block diagram of the basic double-loop sys- 
tem investigated is shown on figure l(a) . The engine, sensors, and con- 
trollers have transfer functions symbolized by E, H, and G. 
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Significant system transfer functions have been derived and are 
presented herein with the added substitution that the product of i 
terms in each simple loop are characterized by one symbol. The product 
of speed loop terms, G^, and is replaced by Ljj, while the 

product of terms in the temperature loop. Eg, Gg, and E^, is replaced by 
Lj. A third loop is formed in this configuration that includes the in- 
teracting engine terms and therefore is called the interaction loop. 

This loop consists of H^, G^, E^, Eg, Gg, and E^. The product of a1 l 

these terms is indicated as in subsequent discussion. 


System transfer functions are 

^ ^[LjjCl+Lp) - Lx] 

Ng (1+Lpj) (l+Lj) - 

J_ ^ ^ 2^5 

Tg (l+Lfj) (1+Lp) - Ly^ 

Tg (1+Lfj) (l+Lj) - 

T Gi®2 

Ng (l+Lp) “ Lx 


( 1 ) 


( 2 ) 


( 3 ) 

(4) 


The stability of the system can be determined from analyis of the 
denominator of these transfer functions, which when set equal to zero 
is the characteristic equation of the system. Further examination of 
this equation, however, indicates that if the interaction loop term Lx 

were made zero, then the system would behave as two independent single- 
loop systems. 

^^^4nteracting system . - A completely noninteracting system can be 
derived by adding two new elements to the control configuration as shown 
on figure l(b). The purpose of X is to add a function of temperature 
error to a function of speed error such that the resulting change in 
fuel flow compensates for the speed change resulting from the action of 
temperature error on exhaust-nozzle area. It follows, therefore, that 
with a properly chosen value of X, no speed error will be evident when 
a change in controlled engine temperature is requested by manipulation 
of set temperature . Another element Y can be added to the system in a 
similar manner so that speed error will have no effect on temperature 
when set speed is varied* 
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Transfer functions for the system shown in figure l(b) have been 
derived and are presented herein. 


N 

^[%.(i+i-rt) - 

%'] 


■ (i4-r^j(i4-i^,) 


N 


3 

T ' 
•^s 

(1-fLjjt ) (1+Lpt ) 

- L^ , 

T 

^CLt' (!+%') - 

%' 3 

T ‘ 
•^s 

(1+Ljji.) (1+Lpt ) 

■ Ix’ 

T 

^1^2 YG2E 

4 

Ns ^ 

(l+Ljjt ) (l+Lpi ) 

- Lj{t 


( 5 ) 


( 6 ) 


( 7 ) 

( 8 ) 


In these equations Lj^i is equal to + YG2E3), which is the 

product of all terms in the speed loop where now a parallel feed path 
exists through and YGgE^. Similarly, L^, is equal to 

H„(G„E. + XG-E-) with the parallel feed being through GgE^ and XG^Eg. 

The interaction loop is given by Lj^t, which is H3_H2(G2_E2 + YG2E4) 

(G2E3 + XG2^E]_) . Two parallel feed paths are evident in this loop. 


The interaction loop can be made equal to zero if either X 

or Y has the following values: 


« 2^3 




( 9 ) 


Y 


QjE2 


( 10 ) 


Substituting equations (9) and (lO) into the characteristic equation for 
the noninteracting system results in the following expression: 


) (1+1*1 ' ) 
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This indicates that the mviltiple-loop system can he considered as con- 
sisting of 2 independent loops. The characteristic equations of these 
loops are 


and the only conditions under which system instability can occur is 
when one of the loops is imstable by itself. 

Further examination of the characteristic equation for the com- 
pensated system shows that only one compensating term X or Y is 
necessary to make the interaction loop eq-ual to zero. However, the 
transient responses of a system with only one added element will be 
different from those of a system with both X and Y added. As an 
example of this, consider only an X element added to the system. 

Speed will not be affected by a disturbance placed in set temperatvire 
even thovigh the system will act to minimize temperature error by caus- 
ing exhaust -nozzle area to change. A disturbance in set speed, however, 
will cause temperature to deviate from initial value as well as cause 
speed to change so as to minimize speed error. The resulting tempera- 
t-ure error, however, will not cause an additional change in speed because 
of the influence of the X term. With the compensating Y element in 
place, tenperature would not be affected by change in set speed. 

An engine control system may not require the complexity of com- 
plete compensation for both temperature and speed error interaction. 
Compensation for the effect of temperature error on speed shotild be suf- 
ficient because normal engine operation is usually at top speed, where 
speed is held constant, and thrust variations are made by requesting 
changes in temperature only. If the compensating element is exactly as 
specified by equation (9), the compensation is complete and no speed 
disturbance results dtiring a transient from a set temperature change. 
However, because speed variation within certain limits can be tolerated, 
it follows that the compensation element need not be so complex as in- 
dicated by eqxaation (9). The analysis reported herein is based on the 
use of only a gain term for the compensating element instead of one 
having all the necessary dynamic terms as indicated by equation (9). 

Figure l(c) is presented to show three possible positions of the 
compensating term in a control configuration. The complete forms of 
X for the three positions can be derived and axe 
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These e:iq)ressions indicate that the con^jensating element will have dif- 
ferent required characteristics depending on the function of temperature 
and speed errors considered. By using only a gain term in the compensat- 
ing element, partial compensation to different degrees is achieved with 

and Xg, while X^ supplies con^lete con^jensation to the system. 

Specific systems investigation . - Figure 2 shows a block diagram 
of the systems investigated as set xxp on the computer using the method 
of reference 3. Component gains and time constants were chosen to be 
representative of current devices, and specific values are shown on the 
figure. 

The engine has a time constant of 1.75 seconds at the operating 
point chosen, which was, based on design values, 96 percent speed, 86 
percent tenperature , 67 percent fiiel flow, and 98 percent exhaust - 
nozzle area, where turbine exit area is defined as 100 percent. Total 
exhaust -nozzle area range is 75 to 133 percent. 

Speed sensor was simulated by a first-order lag having a time con- 
stant of 0.05 second, while tenperature sensor was assumed to be a 
thermoco\aple with a nominal time constant of 1 second. 

The fuel flow servo was represented by two lags in series, each 
having a time constant of 0.10 second. The exhaust -nozzle area servo, 
which in practice is a much slower device, was considered to consist 
of a 0.3-second time constant lag in series with a 0.15-second time 
constant lag. The system was calibrated in such a manner that loop 
gains could be read directly from dial settings. 

fi^st system investigated consisted of the basic configuration 
where speed is controlled by fuel flow and temperature by exhaust- 
nozzle area, with proportional control in both loops. This system was 
then modified by the addition of an integral term to each loop. Integral 
action results in elimination of droop that is characteristic of propor- 
tional controls, The integral time constant was taken to be equal to 
the engine time constant. A third system studied consisted of adding a 
gain term to the basic system to compensate for the effect of tenperatwe 
error on speed. Compensating element Xji^ was used, and as mentioned 

before, this element provides only petrtial compensation. The fourth 
system investigated used both the integral and compensating terms of the 
previous configurations. 

In considering the gross effect of loop gains, it becomes apparent 
that sensitivity of control is related to this gain. A high loop gain 
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results in high sensitivity and rapid recovery to an imposed disturbance. 
The transient responses of this system^ however, become more oscilla- 
tory as loop gain is increased. At some value of this term, depending 
on the dynamics involved, the entire system can become unstable, at 
which point a self* -sustained oscillation will occxir, as shown in refer- 
ence 4. Preliminary investigations were conducted to determine the 
effect on stability limits of the three forms of coii 5 )ensation X^, X^, 

and Xq. 


STABILITY LIMITS 

FigTjre 3 presents stability limits obtained with the basic system 
and also with each of the three gain compensation terms X^^, , and 

Xq. In these data, speed loop gain Ku has the same significance 
whether or not the compensation term is used. When compensation is 
used the gain of the ten 5 >erature loop is actually that computed from 
Lt’, which is 


( ^ 2^5 

^ “ E1E4 

These data, however, are plotted for comparison purposes on the basis 
of the simple temperatiire loop Lp which has a gain symbolized by Kp. 

These conditions account for the coincidence of the stability 
limits of all systems in the high Kjj and low Kp region where the 

system is predominately a single speed-fuel flow loop and for the di- 
vergence of points in the region of low Kjj and high Kp where the 

system is predominately a single ten5)erature-area loop. 

Examination of the curve showing the stability limit with no com- 
pensation (curve labeled K) indicates that the interaction loop has a 
severe effect on the stability of the system in the high temperature 
loop gain region and acts to reduce the gain that this loop could toler- 
ate were it operating alone. The stability limit for the system with 
Xq, in which case only a gain term provides conplete coii 5 )ensation, shows 

that when the product of the terms of the interaction loop is zero the 
stability limit approaches the theoretical limit. The resiat is that 
each loop is independent of the other to the single-loop stability 
limit. The slight deviation of the limit derived using the analog from 
the theoretical limit can be attributed to minor, inaccuracies in ad- 
justment of the compensating gain and too small dynamic terms associated 
with the computer elements. 
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The system en 5 )loying is shown to he more effective than that 

using Xg and therefore was used as the compensation form for subsequent 
work. Xg, in fact, decreased the stability limit below that with no 
compensation. The variations in effect of Xa and Xb can be explained 
if the necessary forms of these crossovers for dynamic compensation are. 
considered. With gain compensation only, X^ is much closer to com- 
plete dynamic condensation than is Xg. Condensation of the form Xg, 
while considered better than X^ or Xg, is unfortunately not useful 

on a real engine system because it is iinpractical to vary fuel flow as 
a f-unction of exha\ist-nozzle area without introducing additional dynam- 
ics to the system. 


Figure 4 shows the stability limits obtained for the four configura- 
tions investigated. These data show that when coinpared with a simple 
proportional control system (curve k) , the addition of the integral term 

^curve K ^1 4 conpresses the stability . limit over the entire region. 

T^e addition of compensation to the proportional control system (curve 
K 4 Kjf) expands the limit in the region of high tenperat-ure loop gains. 


The addition of an integral term to the system with condensation (curve 



compresses the stability limit to a small extent. 


but it 


can be noted that a significant inprovement is still evident when com- 
pared with the limit ctjrve for the proportional -plus-integral 
configuration. 


TRANSIENT STUDIES 

Knowledge of stability limits is not sufficient to characterize 
a system from all points of view. It is necessary to determine how a 
system reacts when some disturbance is introduced. Points of interest 
are related to engine safety, speed of response, and nature of error in 
all pertinent engine parameters dviring transient operation. Transient 
characteristics of the four systems having stability limits presented on 
figure 4 were therefore investigated. In all cases data were obtained 
with a step disturbance placed in set temperature. This disturbance was 
considered to be a 1 percent change in required temperatiire based on 
sea-level rated temperatxire (absolute) as being 100 percent. Data were 
taken at numerous operating points in the stable region of each system. 
Maximum speed and tenperature excursion were recorded. Engine safety 
as related to overspeed and overtemperature can be determined from an 
examination of maxiimum excursion data. 
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A general indication of the speed of responses can be obtained by- 
measuring the time in seconds for the engine to reach maximm speed 
excursion after start of transient, and this time was recorded. 

For the proportional control system, speed and temperat-ure droops 
were also noted. 

These data for the various systems are plotted as contour lines on 
their respective stability limit maps. This presentation permits evalu- 
ation of the effect that either loop has on the other and also enables 
a comparison to be made between systems investigated. Contour lines 
are not extended to the stability limit line because the systems become 
too oscillatory and critical to adjustment in the region close to the 
limit. In addition to these data, photographs of transient responses 
of act\ial temperature Tg^, measured tenrperat\jre speed N, fuel flow 

Wj, and area A were taken at a nvmiber of operating points. On the 

photographs of the transients, amplitude sensitivity of set temperature 
disturbance represents 1 percent change. The same amplitude sensitivity 
applies to all traces. 

Proportional and proport ional-plus- integral controls . - Data pre- 
sented on fig-ure 5 show that for both the proportional and proportional-, 
plus -integral' controls maximum speed excursion is a fxmction of speed 
and temperature loop gains. In both systems, speed excursion decreases 
as speed loop gain is increased and increases as ten^erature loop gain 
is increased. These facts can be explained by the following considera- 
tions: High-speed loop gains resiilt in a sensitive control so that 

small off -speed signals during a transient result in large correcting 
signals which tend to decrease the speed overshoot. However, with in- 
creasing values of temperatiire loop gain, the gain or sensitivity of 
the interaction loop also increases. Therefore a small temperat\are 
error signal diiring the transient introduces a large opposing signal 
into the speed loop , which results in a corresponding Increase in speed 
excursion. 

The system with integral added produces a slightly greater speed 
overshoot d-uring the transient at con^jarable operating points than does 
the proportional control. However, the advantage of this system is that 
no steady-state error or droop exists regardless of loop gains. 

The magnitude of change in droop in the proportional system is 
shown on figure 6. These values were calculated from a consideration 
of equation (l) and were also derived by analog methods. Droop follows 
the same trends as does speed exc-ursion in that it decreases with speed 
loop gain but increases with temperature loop gain. 

Figure 7 presents a conq)arison of the time required to reach the 
point of maximvim speed excursion for both systems. Contour lines of 
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constant times on the stability limit map indicate that at low values 
of tenrperatiure loop gain the times axe very nearly equal. As this loop 
gain is increased the divergence also increases with the proportional 
control being a little faster for the greater part of the range of 
speed loop gain. At high speed loop gains the data show a slight ad- 
vantage to the system with added integral term. This, however, is in 
an undesirable region of control operation because the operating point 
is too close to the stability limit and the system is very oscillatory. 

Maximum excursion of turbine discharge tenperature is presented on . 
figure 8 for the basic configuration and on figure 9 for the system 
with integral added. Temperature data were recorded at two locations 
in the temperature loop. One signal represents actual gas temperature 
T while the other is the thermocouple output or measured temperature 

£L 

Tjn. Under practical conditions the thermoco\:iple indication is the more 

realistic one to use because it is the actual control parameter and also 
it offers a better indication of turbine blade temperature. When operat- 
ing a control system with low loop gains the entire system response is 
slow and a condition of no overshoot or at least very small overshoot 
beyond final value can be established. Under these conditions a thermo- 
couple can follow actual gas temperature with reasonable accuracy. 
However^ at higher loop gains this is not true, and a greater divergence 
between actual and measured maximiim temperature excursion can be ex- 
pected. These conditions are shovm on figures 8 and 9. 

Without the integral term it is possible at low values of ten^er- 
ature loop gain to have maximum excursions of temperature that are less 
than the required change. In addition, the final value is always less 
than the required change because of the characteristic droop associated 
with purely proportional control systems. The addition of the integral 
term results in zero steady-state error and the maximum excursion will 
be at least equal to the required change in set temperature, even at 
low values of temperature loop gain. This basic difference of the two 
systems is shown in figures 8 and 9, and these figures also show that 
at comparable operating points of speed and temperature loop gains the 
maximum exciirsion is greater for the system having the integral terms 
included in the configuration. 

Temperature droop for the proportional system is shown on figure 
10 and indicates that droop decreases as temperature loop gain increases. 
Increases in speed loop gain also tend to decrease droop, but to a 
lesser degree. 

Photographs of significant traces are shown on figure 11 for the 
proportional control system and on figure 12 for the proportional-plus - 
integral system. Examination of these photographs in conjunction with 
data already presented will indicate the magnitude and nature of 
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transient responses at various operating points of the systems. Figure 
ll(c) taken with a speed loop gain Kjj of 1 and a temperature loop gain 

Kp of 0.5 shows that the proportional system is very stable with small 

overshoots^ but that it is inherently slow in response and has a droop 
in both speed and temperature. Figure 11(d) taken with Kt increased 
to 2 shows that the system now becomes more oscillatory with a relatively 
low frequency of superimposed oscillation. Temperature droop is notice- 
ably reduced. Figure 11(a) presents the conditions when is again 

set at 0.5, but % increased to 10. These responses indicate a much 
faster system then present in figure 11(c), but fuel flow and actual 
temperature excursions are greater. The increased actual temperature 
overshoot, however, is of such short duration that it does not contribute 
significantly to the maximum excursion of measured temperature, which 
is more nearly representative of the manner in which turbine blades 
respond. 

Figure 11(b) presents responses taken with Kjj set at 18 and Kp 

at 2. These responses indicate two modes of oscillation before stable 
operation is achieved. Investigation of this action shows that the 
lower frequency is due primarily to the temperature loop which contains 
the slower servos, while the higher frequency is due to action of the 
speed loop which includes faster servos. Actual values of superimposed 
frequencies are not directly determinable from consideration of each 
loop independently because of the effect of the interaction loop. 

A comparison of figure 12 with figure 11 shows that the addition 
of the integral term does not alter the trends indicated by the propor- 
tional system. The two points of difference are that with the integral, 
the droop in speed and temperature is eliminated and the over-all system 
responses axe slower and more oscillatory. 

Controls with added compensation for interaction . - The investiga- 
tion was continued with an analysis of the transient response charac- 
teristic of proportional and proportional -plus -integral systems after a 
compensation term was added, as shown on figure 2. Data indicated 

that these two compensated systems followed similar trends in regard to 
the characteristics of responses; therefore subsequent discussion will 
be based on the compensated Integral system. The only significant dif- 
ference is that the compensated proportional system has a temperature 
droop which is predictable from consideration of equation (7). No speed 
droop is obtained when a disturbance is placed in set temperature because 
the compensation term is so designed that no steady-state speed change 
will result from that disturbance. The system will have a speed droop 
if the disturbance is placed elsewhere in the configuration. 
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The compensated proportional -plus-integral system is a little 
slower in response than the one without the integral term^ hut the ad- 
vantages of the integral action in eliminating steady-state speed and 
temperature errors regardless of where disturbance occurs makes the 
integral action more attractive. 

Figure 13 shows maximum speed excursion data and time to reach 
this peak point for a disturbance in set temperature. Maximum speed 
excursion increases with increasing temperature loop gain and decreases 
with increasing values of speed loop gain. However, it can be noted by 
comparison with figure 5, which shows the function for a noncompensated 
system, that the addition of the con5>ensation term greatly reduces the 
speed loop gain effect on the system, and in addition that the magnitude 
of peak error is greatly reduced at comparable loop gain points. This 
indicates the effect of the compensation term in the system. The small 
speed loop gain effect would be eliminated completely if the compensa-- 
tion term had incorporated in it the necessary dynamic characteristics 
as required by equation (9) . 

Contour lines of time to reach maximum speed excursion point appear 
to follow the general shape of the stability limit. Comparison of these 
data with figure 7 shows that the compensated system is much faster than 
the noncompensated control. 

Turbine discharge temperature characteristics are presented on 
figure 14. With low temperature loop gains and over the full range of 
speed loop gains, no overshoot in actual temperature occurs, so that 
the maximum temperature excursion becomes equal to the' required value. 
This temperature change to final value is primarily due to integral 
action in the system. The same effect can be observed in measured tem- 
perature data but continues to higher values of temperature loop gain 
because of the inability of the thermocouples to follow overshoots in 
temperature. At higher temperature loop gain values the data show that 
temperature excursion is dependent on and increases with teirperature' 
loop gain. 

These data also show that actual and measured temperatures are 
practically. independent of speed loop gain up to the limiting lines 
shown on the map and designated "limit Tg^" and "limit Tj^." At speed 

loop gains above this limit a pronounced dependency does exist. Exami- 
nation of figure 15, which consists of -photographs of typical transient 
responses, will serve to define the nature of these limits. Maximum ex- 
cursion of actual temperature on figure 15(d)' occurs on the first peak 
of the oscillatory response; while on figure 15(b) it occurs on the 
second peak. The values of loop gains which result in equal anplitude 
of the first and second peaks define the limit shown. 
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Comparison of responses shown on figure 15 with those of figure 12 
show that compensation produces very desirable improvements in system 
performance in the region of operation defined by figures 15(c) and 15(d) 
because the system responds faster with less superin^^osed oscillation. 
Further examination of figure 15 shows that two modes of oscillation 
occur at the higher speed loop gain settings as shown on figures 15(a) 
and 15(b) . The higher frequency is the result of interaction of the 
speed loop with temperature loop. It follows therefore that when speed- 
loop gain is low the assumption . can be. made that the simple gain com- 
pensation for interaction is sufficient to allow analysis based on 
single loop considerations. However^ at high values of speed loop gain 
this assumption is no longer valid and additional compensation for dy- 
namic terms is required if it is desired to make the two basic loops 
independent of each other. 


CONCLUDING REMARKS 

Addition of integral terms to the speed and temperature loops com- 
presses the stability limits and makes the system slower and more 
oscillatory than would be the case with proportional control only. How- 
ever., integral a.ction^ by eliminating droop ^ becomes desirable for control 
application because it makes final values of system parameters independent 
of the source of disturbance. 

Addition of complete compensation (dynamic compensating terms) for 
engine interaction eliminates the effect of system interaction^ which 
then permits each loop to be developed individually for a desired re- 
sponse. This compensation has the particular advantage that the compli- 
cated double-loop interacting system has been reduced to two noninteract- 
ing single loops^ and the analysis and synthesis procedures of a single- 
loop servo theory can be applied. 

Addition of proper partial compensation (gain compensating terms) 
results in considerable improvement in the characteristics of an inter- 
acting control system^ and in a practical sense^ is considerably easier 
to apply to a system than complete compensation. 

In particular^ it was found that with the engine operating near 
top speedy where an increase in thrust is obtained by increasing the 
temperature; a single partial compensation tenn from tenperature error 
to speed error resulted in appreciable improvements in system character- 
istics. The system was more stable^ and also faster response times were 
observed. These improvements in characteristics can be considered to 
advantage in a coniparison with the noncompensated system. At comparable 
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gain points, the compensated system is not only faster, but has a 
definitely larger margin of gain to instability. It also follows that 
for comparable responses of the two systems, the requirements on re- 
sponse of the control servos need not be so severe when compensation 
is used. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics • 
Cleveland, Ohio, December 21, 1953 
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APPENDIX - SYMBOLS 


General Symbols 


area of variable-area exhaust nozzle 


gain of fuel flow to speed control loop 


engine gain of speed to area 


engine gain of speed to fuel flow 

gain of area to temperature control loop 

engine gain of turbine -discharge ten^ierature to area 

engine gain of turbine -discharge teniperatvire to speed 

engine gain of turbine-discharge ten 5 )erature to fuel flow 

gain of compensation term 

actual engine speed 
speed error, Ng - H^N 

meastired engine speed, H^N 
desired engine speed 
con 5 )lex Laplacian operator 

actml turbine-discharge tenqperature (Tg^ used when differentiating 
from Tjj^) 

tengjeratvire error , Tg - H 2 T ' 

measured turbine -discharge ten 5 >erat'ure , ' 

desired tixrbine-discharge temperature 


■ engine fuel flow 
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Trtosfer Functions 

Ej_ speed to fuel flow 

Eg terr^Derature to fuel flow 
E^ speed to area 

E^ temperature to area 

Gj fuel, flow controller 

Gg area controller 

speed sensor 

Hg temperature sensor 

X complete condensation term from temperature error to speed error 
partial compensation term from temperature error to speed error 

Xg partial compensation term from area to speed error 

X^ partial condensation term from area to fuel flow 

Y complete con 5 >ensation term from speed error to temperature error 
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(c) Three forms of partial compensation. 
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Figure 2. - Block diagram of system on analog computer. 
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Temperature loop gain, Kji 


Figure 5. - Maximum speed excursion characteristics of proportional control 
compared with proportional- plus- integral control. 




Speed droop, percent (N error is negative) 
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Figure 6. - Speed droop in proportional control with dis- 
turbance in set temperature. 
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Figure 7. - Time to reach maximum speed excursion with proportional control 
compared with time when using prpportional-plus- integral control. 
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(c) 0.5; 1.0. (d) Krp, 2; 1. 

Figure 11. - Proportional control. Transient response to disturbance in set temperature. 
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(a) 0,5; Kjj, 10. 
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(d) 1.0; 1.0. 


Figure 12. - Proportional- pi us -Integral control. Transient response to disturMnce 
in set temperature. 
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Figure 13. - Maximum speed excursion and time to maximum speed excursion for noninteracting 
proportional-plus- integral control superimpos-ed upon stability map. 
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Figure 14. - Actual and measured maximum temperature excursions for non-interacting 
proportional- plus -integral control superimposed upon stability map. 
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Figure 15. - Nonlnteractlng proport Iona 1-plua- integral control. Transient response to 
disturbance In set temperature. 
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